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The photo-nuclear production of vector mesons in ultraperipheral heavy ion collisions is investi- 
gated within the collinear approach using different parameterizations for the nuclear gluon distribu- 
tion. The integrated cross section and the rapidity distribution for the AA — ■> V AA (V = J/^, T) 
process are computed for energies of RHIC and LHC. A comparison with the recent PHENIX data 
on coherent production of J/^f mesons is also presented. We demonstrate that the study of the 
exclusive quarkonium photoproduction can be used to constrain the nuclear effects in the gluon 
distribution. 

00 . 

PACS numbers: 12.38.-t; 25.75.-q; 13.60.Le 

o ■ 

<n : 



43 



A systematic measurement of the nuclear gluon distribution is of fundamental interest in understanding the parton 
structure of nuclei and to determine the initial conditions of the quark gluon plasma (QGP) predicted to be formed 
in central heavy ion collisions (See e.g. [![). Another important motivation for the determination of the nuclear gluon 
^| ' distribution is that the high density effects expected to occur in the high energy limit of QCD should be manifest in the 
modification of the gluon dynamics 0] ■ At the moment the behavior of this distribution is completely undetermined 
by the fixed target experiments, with a possible improvement in the future electron- nucleus colliders (See e.g. [3,0]). 
However, as the date of construction and start of operation of these colliders is still in debate, we need to obtain 
I 1 alternative searches to estimate the medium effects in the nuclear gluon distribution. It has motivated several authors 
Qh) to propose the study of different observables in distinct processes to constrain the nuclear gluon distribution (For 
a recent review see Ref. Q). In particular, in Ref. 6] it was proposed to study the vector meson production in 
ultraperipheral heavy ion collisions at RHIC and LHC in order to constrain the nuclear medium effects present in 
the nuclear gluon distribution. The basic idea is that in these collisions the hi gh flux of quasi-real photons from 
one of the nucleus provides a copious source of photoproduced reactions 0, H, FqI. Il0j|. A photon stemming from the 
electromagnetic field of one of the two colliding nuclei can penetrate into the other nucleus and interact with one or 
more of its hadrons, giving rise to photon-nucleus collisions to an energy region hitherto unexplored experimentally. 
For example, the interaction of quasi-real photons with protons has been studied extensively at the electron-proton 
collider at HERA (For a recent review see Ref. [Hj])- The obtained 7p center of mass energies extends up to 
W-yp ~ 300 GcV, an order of magnitude larger than those reached by fixed target experiments. Due to the larger 
number of photons coming from one of the colliding nuclei in heavy ion collisions, a similar and more detailed study 
qq ' will be possible in these collisions, with W 1 n reaching 950 GeV for the Large Hadron Collider (LHC) operating in its 
heavy ion mode 

ft! HUD- As the cross section for the diffractive vector meson production depends (quadratically) 
on the gluon distribution, it gives a unique opportunity to study the low x behavior of the gluons inside the nucleus. 
The results from Ref. [|| demonstrate that the study of photoproduction of vector mesons in ultraperipheral collisions 
determines the behavior of the nuclear gluon distribution in the full kinematic region, with the rapidity distribution 
allowing for the first time to estimate the EMC, antishadowing and shadowing effects. However, these results were 
obtained considering several approximations for the calculations of the diffractive vector meson production, the GRV95 
parameterization [l2| for the nucleon parton distribution and the EKS parameterization (l3j for the nuclear effects. 
All these models have been improved recently. Furthermore, the PHENIX Collaboration recently released the first 
(preliminary) data on the cross section of the coherent J/ 1 ^ production in gold - gold ultraperipheral collisions at 
y/s = 200 GeV [Til ], providing the first opportunity to check the basic features and main approximations of the 
distinct approaches describing nuclear vector meson photoproduction. Finally, in near future the Large Hadron 
Collider (LHC) at CERN will start its experimental physics program and coherent hadron-hadron interactions should 
be studied [15| . These facts motivate a revision of the estimates presented in Ref. In particular, in this paper we 
revise the predictions for J/ty production in coherent AA collisions and present, for the first time, predictions for the 
T production in coherent collisions considering the collinear formalism. 

Over the past few years a comprehensive analysis of vector meson production in ultraperipheral heavy ion collisions 
was made considering different theoretical approaches [a, lis llZI, 1 1 8L 1 1 9l. 1 2 OL 1 2 ll . |22| . In particular, much effort has been 
devoted to obtain signatures of the QCD Pomeron in such process [12, lU, [23| , which can be used to constrain the QCD 
dynamics at high energies. In this paper we restrict our study of vector meson production to the collinear formalism, 
based on the QCD factorization theorem [24| , which predicts that at small- 2 and for a sufficiently hard scale, the cross 
section of coherent production of vector mesons off any hadronic target, including a nucleus, is proportional to the 
square of the gluon parton density of the target. To lowest order, the 7/1 — > Vh (h — p, A) amplitude can be factorized 
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FIG. 1: (Color online) Energy dependence of the exclusive J/ty (left panel) and T (right panel) photoproduction cross section. 
Comparison with HERA data HE, [H, [33] . 



into the product of the 7 — > QQ transition (Q = c, b), the scattering of the QQ system on the target via (colorless) 

two-gluon exchange, and finally the formation of the quarkonium from the outgoing QQ pair. The heavy meson mass 

My ensures that perturbative QCD can be applied to photoproduction. The calculation was performed some years 

2 

ago to leading logarithmic (log(Q )) approximation, assuming the produced vector-meson quarkonium system to be 
nonrelativistic (25l . |26| and improved in distinct aspects (27l . |28| . Assuming a non-relativistic wave-function for the 
vector meson one have that the t = differential cross section of photoproduction of heavy vector mesons in leading 
order collinear approximation is given by [25L [26j 

t=o = — "73 — - _s X [x9h{x,Q ) , (1) 

at 48a q 

where xgh is the target gluon distribution, x = AQ /W 2 with W the center of mass energy and Q = My /A. Moreover, 
r ee is the leptonic decay width of the vector meson. The strong dependence on xg offers an opportunity to use the 
experimental HERA data for the 7*7? — > J/^p process to determine the behavior of the gluon distribution at low x 
and in the Q 2 region, which is not constrained by the global analyzes [2i|. In Refs. [27l. l28l [29l. l30l. l3l| the authors 
have estimated the relativistic corrections [0(4%)] , the real part contribution of the production amplitude [0(15%)], 
the effect of off-diagonal partons [0(20%)] and next-to-leading order corrections [0(40%)] to the LO exclusive heavy 
vector meson production, given by Eq. (fT]). 

In order to obtain a baseline for our calculations of the vector meson production in interactions we initially 
estimate the 7p cross section and compare it with the HERA data. Following Refs. [l^,[3l| we estimate the total cross 
section for the J/^f and T photoproduction at HERA using the MRST(LO) parameterization [32j for the nucleon 
gluon distribution and including the corrections discussed above. Assuming an exponential parameterization for the 
small |t| behavior of the amplitude, one have that a (jp — > Vp) — \ da (i h ^* vh ) \ t _ Qj where b is the experimentally 
measured slope parameter. In our calculations we assume b = 4.5 GeV -2 as in Refs. [2?l [2^|. which is in agreement 
with [331 ]. For simplicity, we assume the same value of b for J '/vp and T production. It is a reasonable approximation, 
since the HERA data shown that b decreases to a universal value of about 4-5 GeV -2 as the scale Q 2 + m v increases, 
independently of the vector meson studied (For a recent review see, e.g., [34|). In Fig. [T]we compare our results with 
the J/^P (left panel) and T (right panel) HERA data [35j, [36|, [37| and demonstrate that it are reasonably described. In 
the case of T production, our prediction is systematically somewhat below the scarce experimental data. In principle, 
the agreement could be improved by choosing a smaller value of b. In what follows we will assume that the corrections 
for the LO exclusive heavy vector meson production are independent of the target, i.e. we consider that it are the 
same in the proton and nuclear case. Moreover, we assume that in the case of nuclear targets, b is dominated by 
the nuclear size, with b ~ R\ and the non-forward differential cross section is dominated by the nuclear form factor, 
which is the Fourier transform of the nuclear density profile. Therefore, we consider that the total cross section for 
vector meson production in 7^4 interactions is given by 

^V^. ^-^ uf.W, (2) 
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FIG. 2: (Color online) Ratio R g = xg A /A.xg N predicted by the DS HH, EKS HKN [U and EPS El parameterizations 
at Q 2 = 2.5 GeV 2 and A = 208. 



(My/2uj) 2 and -F(t) = f d 3 r p(r) exp(iq • r) is the nuclear form factor for the distribution. Here we 
use the analytical approximation of the Woods-Saxon distribution as a hard sphere,with radius Ra, convoluted with 
a Yukawa potential with range a = 0.7 feu. Thus, the nuclear form factor reads as [161 ]. 



F{q = VW\) = -rr [sHqRa) ~ qR A cos{qR A )} 
Aq^ 



1 



1 + a 2 q 2 



(3) 



where pq = 0.16 fm~ 3 . 

The main input in the calculations of the quarkonium production cross section in photon-nucleus interactions 
is the nuclear gluon distribution. In the last years several groups has proposed parameterizations for the nuclear 
parton distributions (nPDF), which are based on different assumptions and techniques to perform a global fit of 
different sets of data using the DGLAP evolution equations [H, [H, HI, SO, S3- Lets present a brief review of these 
parameterizations (For details see, e.g., Ref. Q). Initially, let's consider the EKS parameterization [l3| which is based 
on the DGLAP evolution equations [42J with adjusted initial parton distributions to describe experimental results 
[43l Hi| presenting nuclear shadowing effects. The basic idea of this framework is the same as in the global analyzes 
of parton distributions in the free proton: they determine the nuclear parton densities at a wide range of x and 
Q 2 > Qo — 2.25 GeV 2 through their perturbative DGLAP evolution by using the available experimental data from IA 
DIS and Drell-Yan (DY) measurements in pA collisions as constraint. In this approach, the nuclear effects are taken 
into account in the initial parton distribution xf A (x, Qq) of the DGLAP evolution. EKS have expressed the results in 
terms of the nuclear ratios Rf(x, Q 2 ) for each parton flavor / in a nucleus with A nucleons (A > 2), at 10~ 6 < x < 1 
and 2.25 GeV 2 < Q 2 < !0 4 GeV 2 . Other groups have considered different sets of data, next-to-leading order (NLO) 
corrections to the DGLAP equation and uncertainty estimates, and have proposed distinct approaches to describe 
the nuclear effects [HI, Htl H3, IH| • Next-to-leading order (NLO) corrections to the DGLAP evolution are taken into 
account in (38l. |40|. In particular, De Florian and Sassot (DS) 38'] proposed a framework where each nuclear parton 
distribution is described by a convolution of the corresponding free nucleon parton distribution with a simple flavor 
dependent weight function that takes into account the nuclear effects. On the other hand, the authors of [39l l40| have 
calculated the nPDF uncertainties by using the Hessian method and have obtained LO and NLO |40l ] nuclear 
parameterizations. More recently, an important step toward better constraining the gluon nuclear distributions was 



done in 41(, where the authors include recent RHIC data from high-pT hadron production at forward rapidities in 
d+Au collisions, which probe small values of x in the gluon distribution of the Au target, in their global analysis. 
The striking result is a much stronger gluon shadowing than in the previous parameterizations. Due to the scarce 
experimental data in the small- a; region and/or for observables strongly dependent on the nuclear gluon distribution, 
the current status is that its behavior is completely undefined. It is demonstrated by the analysis of the Fig. [21 
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FIG. 3: (Color online) Rapidity distribution for J/$ photoproduction in AuAu collisions at RHIC. 



where we present the results for the ratio R g = xg A /A.xg N predicted by the EKS [13], DS [II], HKN [Io| and EPS 
[4~fl ] parameterizations at Q 2 — 2.5 GeV 2 and A = 208. As we can see, these parameterizations predict very distinct 
magnitudes for the nuclear effects. For larger values of x, the EKS and the EPS show antishadowing, while this effect 
is absent for the HKN and EPS parameterizations in the x < 10 _1 domain. The more surprising feature is however 
the amount of shadowing in the different parameterizations. While the shadowing is moderate for DS and HKN 
parameterizations and somewhat bigger for EKS one, the EPS prediction has a much stronger suppression compared 
with the other parameterizations. For smaller x around x ~ 10 -5 , while DS and HKN parameterizations have about 
20% suppression and EKS one have about 40% suppression, for the EPS parameterization this effect goes to almost 
80% suppression in the nuclear gluon compared with the A scaled gluon content in the proton! For bigger values 
of x the behavior is distinct for all parameterizations. As x grows, the DS parameterization predicts that R g grows 
continuously to 1, that means that the shadowing dies out when x — > 10 -1 . The same happens for the HKN one 
in this limit, but this growth starts only at x > 10 -2 , with R g being flat for 10 -5 < x < 10 -2 . At x « 10" 1 , we 
have that behaviors predicted by the EKS and EPS parameterizations are similar, with R g exceeding 1.2. The main 
distinction between these parameterizations is that in the EPS parameterization one has a much steeper growth, from 
a much stronger suppression at smaller x to the antishadowing behavior for larger values of x. The difference between 
the distinct parameterizations observed in Fig. [5] will be amplified in the exclusive quarkonium photoproduction in 
"fA interactions due to the quadratic dependence on xgA of the cross section. Consequently, it is expected a large 
difference between the predictions obtained using the distinct parameterizations. This is the main motivation for our 
calculations. 

We now present a brief review of vector meson production in ultraperipheral relativistic heavy-ion collisions (UPC's), 
introducing the main formulas. In UPC's the ions do not interact directly with each other and move essentially undis- 
turbed along the beam direction. The only possible interaction is due to the long range electromagnetic interaction 
and diffractive processes (For a review see, e. g. Refs. @, H, H, EH)- In the particular case of photonuclear reactions 
the cross sections are given by the convolution between the photon flux from one of the nuclei and the cross section 
for the scattering photon-nuclei. The final expression for the production of vector mesons in ultraperipheral heavy 
ion collisions is then given by, 

oo 

OAA^AAV {^J SWn) = J du> — &>yA->VA (w^ A — 2b)y/ Snn) (4) 

where V — J/^P or T, uj is the photon energy with to m i n — my jA^jjm v and V^nn ^ s ^ ne io n -fo n c.m.s energy. The 
Lorentz factor for LHC is jl = 2930, giving the maximum c.m.s. jN energy W^a ^ 950 GeV. In this process the 
nuclei are not disrupted and the final state consists solely of the two nuclei and the vector meson decay products. 
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TABLE I: The differential cross section at midrapidity (da/dY\ y —o) for nuclear vector meson photoproduction in UPC's at 
RHIC and LHC energies. At RHIC we also present in brackets the predictions considering the nuclear breakup probability. 



Consequently, the final state is characterized by a small number of centrally produced particles, with rapidity gaps 
separating the central final state from both beams. Moreover, due to the coherence requirement, the transverse 
momentum is limited to be small. Therefore, these reactions can be studied experimentally by selecting events with 
low multiplicity and small total px- It is important to emphasize that recent STAR p° photoproduction data (45j 
show significant coherent production to pr < 150 MeV. The photon flux is given by the Weizsacker- Williams method 
Q • The flux from a charge Z nucleus a distance b away is 



d 3 N (u, b 2 ) _ ft_ 



diu d 2 b 7r 2 uib 2 



(5) 



where 7l is the Lorentz boost of a single beam and rj = tob/^/L] Kq^i(j]) are the modified Bessel functions. The 
requirement that photoproduction is not accompanied by hadronic interaction (ultraperipheral collision) can be done 
by restricting the impact parameter b to be larger than twice the nuclear radius, Ra — 1.2 A 1 / 3 fm. Therefore, the 
total photon flux interacting with the target nucleus is given by Eq. ([5]) integrated over the transverse area of the 
target for all impact parameters subject to the constraint that the two nuclei do not interact hadronically. An analytic 
approximation for A A collisions can be obtained using as integration limit b > 2 Ra, producing 



dN(uj) 2Z 2 a em [ rj 2 . 2 2 , 

- — — — V K Q (ry) Ki (77) + — [K{ (77) - K a (77)) 



duj 



(6) 



where fj — 2uj RaHl- It is worth mentioning that the difference between the complete numeric and the analytical 
calculation presented above for the photon flux is less than 15 % for the most of the purposes [71 la. 19. Il0|. 

In what follows we calculate the rapidity distribution and total cross sections for quarkonium production in 
ultraperipheral heavy ion collisions at RHIC and LHC energies. Moreover, since none of the RHIC and LHC 
detectors have full rapidity coverage, we also present our predictions for dcr/dY at Y = 0. We assume that 
xgA(x,Q 2 ) = R„(x,Q 2 ).Axgn(x,Q 2 ), with R g given by the DS, EKS and EPS parameterizations and xg p given 
by the MRST(LO) parameterization [32|. As the HKN parameterization is similar to the DS one in the kinematical 
range considered, it is not included in our analyses. The distribution on rapidity Y of the produced final state can be 
directly computed from Eq. by using its relation with the photon energy cu, i.e. Y oc In (2Lo/my). Explicitly, the 
rapidity distribution is written down as, 

da [A + A-f A®V ® A] dNJu) . , 

— dY ^«-^a 7 A^(«) (7) 

where ® represents the presence of a rapidity gap. Consequently, given the photon flux, the rapidity distribution is 
thus a direct measure of the photoproduction cross section for a given energy. Initially we present our predictions 
for J/\I> production in coherent AA collisions at RHIC, considering A — Au and ^/smn = 200 GeV. The production 
at mid-rapidity at the RHIC probes ^-values of order 10 -2 , where there is a large difference between the distinct 
nuclear parameterizations (See Fig. |2|). For comparison we also present the prediction obtained assuming R g = 1.0 
(no shadowing). Our results are shown in Fig. [31 We have that the DS prediction is similar to the no shadowing 
one. This behavior is expected due to the small magnitude of the nuclear effects present in this parameterization at 
x > 10~ 2 . On the other hand, the presence of the antishadowing effect in the EKS and EPS parameterizations modify 
the rapidity distribution at large \ Y\, implying an enhancement in the distribution. The study of this kinematical range 
can be useful to constrain the magnitude of the antishadowing effects. In contrast, the behavior of the distribution at 
mid-rapidity is directly associated to the magnitude of the shadowing effect in the nuclear gluon distribution, which 
is larger in the EPS parameterization in comparison to the EKS one. It implies the larger suppression of the rapidity 
distribution at Y — when calculated with the EPS gluon observed in Fig. [3j These behaviors are observed in the 
predictions for da/dY \ y=0 show in Table Q] 

Recently, the PHENIX collaboration has released the first (preliminary) data on the differential cross section of 
the coherent J production at mid-rapidity in gold - gold ultraperipheral collisions at \fs — 200 GeV accompanied 




FIG. 4: (Color online) The rapidity distribution for nuclear vector meson photoproduction on UPC's in AA reactions at LHC 
energy (^S^n = 5.5 TeV). 



by Au breakup The value measured was: da/dY\Y=Q — 48. ± 14.(stat) ± 16.(syst) fj,b. In order to compare 

with this PHENIX data we scale down our predictions by the nuclear breakup probability Pxn ~ 0.64, which takes 
into account the reduction of the yield expected when requiring coincident forward neutron emission. It corresponds 
to J production in coincidence with Coulomb breakup of at least one of the nuclei (For a detailed discussion see 
Refs. 0, When compared with the preliminary PHENIX data we have that it is reasonably described by 

the EPS prediction (See Table H}. Although it is a very interesting result, since the typical values of x probed at 
mid-rapidity are similar to those at forward rapidities in central d+Au collisions, which have motivated the inclusion 
of a stronger shadowing in the EPS parameterization, the current uncertainties and low statistics preclude yet any 
detailed conclusion. As shown in Ref. [l5[ several models, based on very distinct assumptions, also describe this data. 
However, differently from the models proposed in Refs. 0, [l?], QjJ [HJ , the EPS result predicts an enhancement of 
the rapidity distribution at large which can be considered a signature of this model. 

In Fig. |4]we present our predictions for J/vI/ (left panel) and T (right panel) production in coherent AA collisions, 
considering A — Pb and ^/snn = 5.5 TeV. The production at mid-rapidity at the LHC probes x-values of order 
(6 — 20) x 10~ 4 , where the nuclear parameterizations differ by a factor 4 (See Fig. [5]). The no shadowing prediction 
is presented for comparison. At LHC energies the rapidity distribution is strongly dependent of the magnitude of the 
shadowing effects in the nuclear gluon distribution. While the DS prediction implies a small reduction at mid-rapidity 
in comparison to the no shadowing one, the rapidity distribution is suppressed by a factor 10 if calculated using 
the EPS gluon. Moreover, the difference between the predictions of the three parameterizations is very large, which 
implies that the rapidity distribution for exclusive quarkonium photoproduction in UPC's at LHC could be used to 
constrain the magnitude of the shadowing effect. 

Finally, in Table Q] we present our predictions for the integrated cross section at RHIC and LHC energies considering 
as input the distinct nuclear gluon distributions. In the RHIC case we present our results considering the suppression 
associated to the probability of nuclear breakup. In this case the no shadowing and DS predictions for the total cross 
section are very similar, as expected from Fig. [31 while the EKS one predicts an enhancement of 10 % and the EPS 
one a suppression of 20 %. On the other hand, at LHC we have that the total cross section is strongly suppressed by 
the shadowing effects. For example, in J (T) production, the EPS prediction is a factor ss 5 (8) smaller than the 
no shadowing one. 

In what follows we compare our results for the integrated cross sections with other theoretical predictions currently 
available. The photoproduction of J/^ in UPC's was estimated in Refs. l2ll [2^ | . In general our predictions 

at RHIC energies are larger than the results presented in those references when the nuclear breakup is included, with 
our EPS prediction being almost identical that shown in Ref. This result is directly associated to the inclusion of 
the antishadowing effects in our predictions and/or due to the different treatment of t he p hoton-nucleus interaction. 
At LHC energies, our predictions are similar those presented in Refs. [l(| [H, H, HJ [Hj]. In particular, our no 
shadowing prediction is almost identical to the impulse approximation prediction pres ented in Table 2 from Ref. [l(| • 
Moreover, our EPS prediction is almost identical to the result obtained in Ref. [17] usin g a leading twist model of 
nuclear shadowing, while our EKS prediction is similar to the results presented in Refs. Il6l. Il9l . w\ . On the other 
hand, the photoproduction of T in UPC's at LHC energies was estimated in Refs. [HI, |20|. l2l|. Our no shadowing 
prediction is almost identical to that presented in Ref. [l8|, where the nuclear effects are disregarded. On the other 
hand, our EPS prediction is smaller by a factor > 3 than the results presented in Refs. [20|, [21j, with the EKS 
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TABLE II: The integrated cross section for nuclear vector meson photoproduction at UPC's at RHIC and LHC energies. 

one being larger by a factor > 1.25 than the predictions from Refs. [IS HI]]. These results indicate that the study of 
T production can be useful to discriminate between the different models. 

As a summary, in this paper we have calculated the rapidity distribution and integrated cross sections of exclusive 
photonuclear production of vector mesons in ultraperiphcral heavy ion collisions within the QCD collinear approach 
in the leading logarithmic approximation. The main characteristic of this approach is the quadratic dependence 
on the gluon distribution, which makes it an excellent probe of the behavior of this distribution. We demonstrate 
that this model reasonably describes the current data on vector meson photoproduction in scattering on protons. 
We extend this model for interactions and estimate the cross sections for the A + A — > A + A + V (V = 
J , T) process for RHIC and LHC energies considering as input different parameterizations for the nuclear gluon 
distribution. We demonstrate that the rapidity distribution and the integrated cross section are strongly dependent on 
the model considered for the nuclear effects. It implies that these quantities are powerful observables to discriminate 
among the different parton distributions in the nuclear medium. In particular, the predicted shadowing by the EPS 
parameterization is considerably larger than in the previous nuclear PDF's. This strong shadowing effect in the 
quarkonium photoproduction could be tested at the forthcoming LHC through the detection of two rapidity gaps 
with the lepton pair from the meson decay identified in the ATLAS, ALICE and CMS detectors. Our main conclusion 
it that the exclusive quarkonium photoproduction in ultraperipheral heavy ion collisions is a useful tool to help in 
obtaining the correct gluon distribution in the nuclear medium. 
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